Theory of enhancement of thermoelectric properties of materials with nanoinclusions 
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Based on the concept of band-bending at metal/semiconductor interfaces as an energy filter for 
electrons, we present a theory for the enhancement of the thermoelectric properties of semiconductor 
materials with metallic nanoinclusions. We show that the Seebeck coefficient can be significantly 
increased due to a strongly energy-dependent electronic scattering time. By including phonon 
scattering, we find that the enhancement of ZT due to electron scattering is important for high 
doping, while at low doping it is primarily due to decrease of the phonon thermal conductivity. 

PACS numbers: 73.63.-b, 72.15.Eb, 72.10.-d, 65.40.-b 



I. INTRODUCTION 

The energy conversion efficiency of thermoelectric de- 
vices depends on the figure of merit ZT = S'^(jT/k, 
where S, a, T, and n are the Seebeck coefficient, elec- 
trical conductivity, temperature, and thermal conductiv- 
ity. In the best thermoelectric materials ZT is typically 
~ 1, and it is difficult to increase ZT beyond this value 
because of competing effects of electrical and thermal 
conductivities. Advances over the past decade show that 
it is possible to enhance ZT in nanostructured thin-film 
systems by taking advantage of quantum confinement to 
enhance the power factor S^a^ or to increase phonon 
scattering at interfaces to reduce the lattice contribution 
to K^ On the other hand, many existing and envisioned 
thermoelectric applications will require a material that 
is itself of macroscopic dimension. Therefore, recent re- 
ports of property enhancement in bulk alloys possessing 
nanometer-scale compositional modulations have gener- 
ated much excitement. "^'^'^'^ ZT values as high as 2.2 have 
been reported^ in the (PbTe)2:(AgSbTe2)i-a; system, 
and have been ascribed to a large Seebeck coefficient and 
low lattice thermal conductivity due to nanoscale clus- 
tering of Ag and Sb. Heremans et sA^ reported that the 
Seebeck coefficient in bulk PbTe can be increased signif- 
icantly by precipitating a fine distribution of Pb nanoin- 
clusions, and suggested heuristically that the increase in 
Seebeck coefficient originates from an energy-filtering ef- 
fect due to a strongly energy-dependent electronic scat- 
tering time. Kim et ali^ observed an enhancement of 
the thermoelectric properties when ErAs nanoparticles 
of 2.4 nm average diameter were embedded in a InGaAs 
matrix, and ascribed the increase to a reduction in the 
phonon thermal conductivity. Given these observations, 
a general understanding of the role of nanoinclusions in 
enhancing the thermoelectric properties of materials is 
needed, in particular to assess the relative importance of 
electronic and phonon scattering. 

In this paper we present a theoretical model and nu- 
merical calculations of the thermoelectric properties of 
bulk semiconductors containing metallic nanoparticles. 
Our model considers scattering of electrons on the band- 
bending at the interfaces between the semiconductor host 
and randomly distributed metallic islands. This causes 



energy-dependent scattering of electrons, leading to an 
energy filterting effect that increases the Seebeck co- 
efficient. This provides an explicit physical model for 
the proposed energy filtering effect i^ By combining this 
model with a model for phonon scattering on the nanoin- 
clusions, we predict significant enhancement of the ZT 
factor. 

We point out that while the role of metallic nanoinclu- 
sions may appear at first to be similar to that of point de- 
fects for which extensive work has been done, the physics 
is actually quite different. Indeed, in addition to the elec- 
tronic scattering, the phonon scattering on nanoinclu- 
sions occurs in a completely different regime than that 
on point defects, as will be discussed in section IVB. 

The central idea in this paper is illustrated in Fig- 
ure 1. There, spherical metallic nanoinclusions with ra- 
dius R and volume fraction x are randomly distributed 
in a bulk semiconductor material. In general, at such 
metal/semiconductor interfaces, charge transfer between 
the metal and the semiconductor leads to band-bending 
away from the interface, characterized by the electro- 
static potential V{r) (Fig. lb). The presence of this 
potential causes energy- dependent scattering of electrons, 
as illustrated in Fig. Ic. The high-energy electrons are 
unaffected by the potential, but the low energy electrons 
can be strongly scattered. Because the Seebeck coeffi- 
cient depends on the energy derivative of the relaxation 
time d In T [E) /dE at the Fermi energy, this type of en- 
ergy filtering is precisely the prescription to increase the 
Seebeck coefficient of thermoelectric materials. 

Our theoretical model is based on the Boltzmann 
Transport Equation (BTE) within the relaxation time 
approximation. We apply the model to a system of n- 
doped PbTe with metallic nanoinclusions because of the 
availability of experimental data and good understand- 
ing of scattering mechanisms in bulk PbTe^'^'^, although 
the theory can be used for any thermoelectric material. 



II. CHARGE AND HEAT TRANSPORT IN 
BULK PBTE 

In this section we will review the expressionsi^i^ for 
the charge and heat transport in bulk PbTe with n-type 
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FIG. 1: (a) Schematic of the semiconductor host with metaUic 
nanoinclusions. Panel (b) shows an example of the calculated 
potential V(r) and the energy diagram for PbTe at T = SOOA', 
n = 2.5 X 10^Vm-^ Vs = -O.lleV, and R = 1.5nm. Panel 
(c) illustrates the concept of energy filtering: low energy elec- 
trons scatter strongly with the potential, but high energy 
electrons are unaffected. The calculated electronic relaxation 
time for the potential of panel (b) is also shown. 



doping. The valence band of PbTe contains four energy 
minima located at the L points. The energy dispersion 
relation near each minima is usually described by the 
Kane model^ 
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where Eg is the direct energy gap of PbTe, h is the Planck 
constant, and k and m* are the electron wavevector and 
effective mass (at minimum energy point fc = 0, -E = 0) 
along the longitudinal (suffix I) and transverse (suffix t) 
directions of the corresponding L point. For n-type PbTe 
the electron concentration is given by 
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where mjj — 4^/'^(m;*m(^)^/'^ is the density of states effec- 
tive mass in which the 4-fold degeneracy is included, kg 
is the Boltzmann constant, fo{z) is the Fermi function 
written in terms of dimensionless variables z = E/{kBT) 
and zp = Ep/{kBT). Ep is the Fermi energy, and the 
function 7(2:) = z + bz'^, where b = ksT/Eg. In the relax- 
ation time approximation the BTE expressions for elec- 
trical conductivity, cr, Seebeck coefficient, S, and electron 
contribution to thermal conductivity, k^, are^ 
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where m* — 3/{l/m'l+2/m^) is the effective conductivity 
mass, and the average is defined as 
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In bulk PbTe at room temperature the dominant con- 
tributions to the total relaxation time, Tiiuik(z), are scat- 
tering by the deformation potential of acoustic and opti- 
cal phonons, and polar scattering by optical phononsi^i^ 
We also take into account scattering on the short-range 
potential of vacancies although it gives a much smaller 
contribution compared to scattering by phonons. Thus, 
the total relaxation time for bulk PbTe is given by 
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The relaxation time due to polar scattering by optical 
phonons reads^i^ 
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where m*^-y = (m^mj'^)^/'^ is the density of states effective 
mass in the single valley, the function ^'{z) = 1 + 26z, eo 
and £00 are the static and high frequency permittivities, 
and 
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The relaxation time due to scattering by the deforma- 
tion potential of acoustic and optical phonons, and also 
due to scattering on the short range potential of vacan- 
cies can be written generally as^i^ 
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where A = bz{l - Km)h'{z), and B = 86z(l + 
bz)Km/{i'y''^(z)), with the suffix m — a for acoustic 



phonons, m = o for optical phonons, and m = v for 
vacancies. The constants tq,^ and Kj^ are defined as 
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Here Ci is a combination of elastic constants, Eac and Eav 
are the acoustic phonon deformation potential coupling 
constants for conduction and valence bands, Eqc and Eov 
are optical phonon deformation potential coupling con- 
stants for conduction and valence bands, Uyc and 11^ 
are coupling constants of the short range potential of va- 
cancies for conduction and valence bands, a is the lattice 
constant, ujq is the frequency of optical phonons, and p 
is the mass density. Ny is the concentration of vacancies 
calculated from the condition that one vacancy gives two 
charge carriers, N^ = n/2. 

The parameters used for calculation of the relaxation 
times in bulk PbTe at T = 300K are taken from Ref. [1]. 
These parameters are shown in Table HI For calculations 
at different temperatures we assumed the values of these 
parameters to be the same as for T = 300K except for 
Eg and m^ which were linearly interpolated and extrap- 
olated using T = 4.6K and T = 300K values^ with Eg 
saturating for T > AOOK^ 



Parameter 


Value 


Parameter 


Value 


E, 


0.315 eV 


hujo 


0.0136 eV 


m*/mo 


0.0453 


a 


6.461 A 


m* /mo 


0.24 


P 


8.24 g/cm 


eo 


400 


Eac 


15 eV 


too 


32.6 


^oc 


26 eV 


Ci 


7.1xlO"'N/m 


Ka.o 


1.5 


U^c 


3xl0~^^erg cm^ 


K^ 


1.5 



TABLE I; Parameters used to calculate the relaxation times 
for bulk PbTe at T=300K.— mo is the free electron mass. 



III. ELECTRON SCATTERING ON 
BAND-BENDING POTENTIAL OF 

NANOINCLUSIONS 

A. Band-bending potential 

In our model we assume that spherical metallic nanoin- 
clusions with radius R and volume fraction x are ran- 
domly distributed in a n-doped PbTe host material. In 



this section we will calculate the contribution to the re- 
laxation time due to scattering of electrons on the band- 
bending potential at the metal-semiconductor interface. 
For a single nanoinclusion, the electrostatic potential 
V{r) can be calculated by solving the Poisson equation 
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The right-hand-side of this expression is simply the spa- 
tially varying charge [see Eq. ^] calculated by assum- 
ing a rigid shift of the electronic bands with the local 
potential V{r). We solve the Poisson equation with the 
boundary conditions V (oo) — and V{R) ~ Vb- {Vb is 
the potential at the semiconductor/metal interface. The 
value of Vb is fixed for a particular metal, and depends 
on the detailed properties of the interface. However, one 
may consider it to be an optimization parameter provided 
that the physics of the metal/semiconductor interface al- 
lows tailoring of Vb by choosing the metal.) We used 
the fourth-order Runge-Kutta and shooting methods in 
order to solve Eq. pS]) . Figure lb shows an example of 
the calculated potential V{r). 

For small values of Vb or for large r (when the potential 
is screened and small) the right-hand-side of Eq. (jisp can 
be linearized with respect to small V: 
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where A is the screening length. The solution of Eq. H 
is 



V{r) ex -e-''/^. 



(17) 



For degenerate electrons the expression for A takes the 
simple form 
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Due to the large value of the dielectric constant in PbTe, 
eo = 400, the screening length for typical doping con- 
centrations is several times larger then the wavelength 
of electrons on the Fermi surface. For example, for 
n = 5x lO^^cm"'^ one can obtain from Eq. P^ A = llnm 
and kpX — 7, where kp is the wavevector of the electron 
on the Fermi surface. This value of kp\ slowly varies with 
doping (for degenerate electrons we can use Eqs ([2]) and 
(fT8| to obtain the dependence on doping kp\ oc -n}/^). 
On the other hand, a large value of A can lead to an over- 
lap of the band-bending between nanoparticles, which 
may change the bulk carrier concentration. We have re- 
stricted our calculations to a parameter range (doping 
and inclusion volume fraction) where such effects are not 
significant. 



B. Relaxation time for scattering on nanoinclusions 

When the electron scattering on nanoinclusions is 
taken into account the total relaxation time r is 



bulk 



(19) 



where the relaxation time for bulk PbTe is given by Eq. 
([7|), and t^ is the relaxation time due to scattering by 
V{r) at randomly distributed metallic inclusions 



Here 



r ^ =: n^vat. 
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is the concentration of inclusions, at is the electronic 
transport scattering cross-section, and v = dpEp is the 
electronic velocity with p the momentum. 

In order to calculate the transport cross-section in 
a system with nonparabolic energy dispersion, we con- 
sider an electron with momentum p and wave function 
^ (r) — Wp(r)e'P''/'' in the periodic field of the unper- 
turbed PbTe crystal of unit volume. Here Up(r) is the 
periodic Bloch amplitude. As mentioned earlier, the non- 
parabolicity of the electron energy dispersion near the 
conduction band minima is usually described by the Kane 
model^ 
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The isotropic energy dispersion in a form of Eq. (j22|) 
with density of state mass rn*^^ is usually assumed in 
the calculation of the relaxation timei^ii^ii^ The transition 
probability for scattering from state ip^. to state ippf per 
unit time due to a perturbation potential V{r) is given 
by the standard formula of perturbation theoryi^ 
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(23) 



Vp^ 



where dv = (i"^p/(27rfi,)'^. The matrix elements are 

/"VvW^W^pW^^i" ~ f e'^P'P'^''/''Vir)d^r, 

(24) 
where we used the fact that collisions with a small mo- 
mentum transfer dominate scattering on the slow varying 
bend-bending potential, therefore the Bloch amplitudes 
entering Eq. (IM)) are rather close to each other and the 
overlap factor is about unity / Up,{r)up{r)d^r «1. 

Applying Eq. ^^ to Eq. (gS]), the calculation of 
dwfi becomes identical to the calculation of the tran- 
sition probability for scattering of a plane wave e^^''"' '^ in 
a model system described by an equation 



(Ep + V)i^ = Ep^ 



(25) 



with unperturbed Hamiltonian Ep, p = —ihdr, and 
perturbation potential V{r). Applying the operator 1 + 
Ep/Eg to Eq. ([25|) one obtains 
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Here we neglect the commutator term E ^[Ep, V{r)]ip 



E-'[Ep,Vir)y. 
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;V{r) < V{r) (27) 



using the fact that rV{r) is a slow varying function and 
kpX > 1. 

Eq. (p6)l has the form of the usual Schrodinger equa- 
tion that can be used for numerical solution of the scat- 
tering problem with a potential 
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C. Calculation of transport cross-section 

The transport cross-section for scattering on the spher- 
ically symmetric potential is given byiS 

7r 

at =27r I \f{e)f{l-cose)smede, (29) 



where f{0) is the scattering amplitude defined from the 
large-r asymptotic of the wave function 
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(30) 



The wave function tp is a solution of the Schrodinger 
equation (|^5|) with potential Up{r) given by Eq. (|28p . 
It can be expressed as a sum of contributions with differ- 
ent angular momentum 1^ 



V' = ^PKcos0)i?feKr), 



(31) 



1=0 



where Pi are the Legendre polynomials and Rki{r) are 
solutions of the radial Schrodinger equation (fc — p/h): 
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dydrRu +ie- ^-H^ - ^Up{r) ] Rui = 0. 
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(32) 



The large-r asymptotic of Rki (r) has the form 

i?fci(r) (xsin(fcr-/7r/2-f (5;)/r, (33) 

where 5i is the phase shift of Rki{r) relative to the 
potential- free solution. The scattering amplitude f{9) 
can be expressed in terms of the phase shifts asi5. 
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The calculation of the transport cross-section ([291) 
with scattering amplitude ([34]) can be performed in 
the same way^" as the usual cross-section [without the 
(1 — coa9) factor in Eq. P5|) ]. In the integrals 



J Pi {cos 9) Pi' {cos 9) {1 — cos 0) sin 6*^6* that appear in the 



evaluation of the transport cross-section ([29]) only terms 
with I' = IJzLl give nonvanishing contributions <W. After 
integration over 9, Eq. (|29p can be expressed in terms of 
the 6i as 



1=1 



(35) 



We used the following numerical procedure to solve the 
Schrodinger equation ([5^ and calculate the phase shifts. 
The spherical Bessel function ji{kr) is a regular solution 
of Eq. p2]) for r < _R, in the region where Up{r) = 0. 
We used the fourth-order Runge-Kutta method to solve 
Eq. ([32]) for r > i? with boundary conditions such that 
the solution Ri{kr) and its derivative match the spheri- 
cal Bessel function at r = i? . At some large r — r^^x 
we assume that the potential Up vanishes and match the 
solution Ri{kr) and its derivative R'i{kr) to a linear com- 
bination of the spherical Bessel function ji (kr) and spher- 
ical Neumann function yi{kr) (which is another solution 
of Eq. §^ for [/p(r) =0): 



Rl{kr)\r=r 



aji{kr) + Pyi{kr) 



(36) 



with a = {Riy[ - R'lyi) / {jiy'i - j[yi) and /? = -{Rij'i - 
R'iJi)/iJiy'i^ JiVi)- Finally the phase shift 6i is given by 



5i = — 'Aicta,n{j3 /a). 



(37) 



In conjunction with equations (|20p and ([35p this pro- 
vide expressions to numerically calculate r^ for a given 
V{r). An example of the numerically calculated r^ is 
given in Fig. Ic for scattering by V{r) of Fig. lb with 
inclusion volume fraction x — 5%. A simple fit gives 
the dependence ti {E) ^ £'i-39 as shown by the dashed 
line in Fig. Ic. This energy dependence of t^ is much 
stronger than that of r(,„;fc , and leads to enhancement of 
the Seebeck coefficient. 

To obtain an analytical description of r^ (E'), we also 
calculated t^ in the Born approximation by using Fermi's 
golden rule [first term in r.h.s. of Eq. ([23]) ] 



1 p^ dp 

7?°^ ~ 27rfi4 Ze 



\Vp,p\^{l-cos9)sm9d9. (38) 



Here 9 is the angle between initial and final momenta 
p and p'. The expression for rf""'' can be simplified 



by taking the angle integration in Eq. ([38|) for Vpip 
and making the substitutions of integration variables 
t — 2fci?sin I and y — r/R. Finally one obtains 



,B..„(^),^3/.(l + W^-giv^^ (39) 



1 + 2E/Eg X 3a{E, R) ' 



where 
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Numerical tests show that for IVbI ^ CleF the Born 



approximation is valid, Ti{E) 



^Born 



{E), while for 



\Vb\ > O.lel/, Tf°'"^{E) begins to deviate from t,{E) 
calculated from the exact solution of Schrodinger's equa- 
tion ([5^ . Nevertheless, Eq. ([5^ allows us to analyze the 
energy dependence of the relaxation time that is difficult 
to do by using the exact formulas ([20]) and ([35]) . For 
energies E > O.leF, the integral over variable t in (HP]) 
weakly depends on the upper limit of the integration, 
and the function a{E, R) varies slowly with both E and 
R. Thus, we have rf""^" (E) ~ E'^/^, in good agreement 
with the full numerical calculations which yielded a de- 
pendence Ei-^s ^ Comparing the result rf"^" {E) ^ E^/^ 
with the expression ([^0]) (and using v{E) ~ VE) we find 
that the electronic scattering cross-section of the band- 
bending potential depends on energy as E~'^\ this strong 
energy dependence is responsible for the superlinear en- 
ergy dependence of r^ {E) . 



IV. RESULTS AND DISCUSSION 

A. Enhancement of the Seebeck coefficient and 
power factor 

The calculation of the total relaxation time allows us 
to obtain S, a and Kg using the expressions ([3][5]). We 
first consider a specific case by adopting a simple model 
for the interface potential Vb = ^m ^ X + Ef with 
^m the metal workfunction and x the electron affinity, 
and choose ^m — X = — 0.35el/ corresponding to Pb 
nanoinclusions (work function $,„ = A.25eV—) and an 
electron affinity for PbTe XpbTe — 4:.6eV—. Figure 2 
shows the calculated room temperature Seebeck coeffi- 
cient as a function of the doping n and fixed nanoin- 
clusion volume fraction x = 5%. We note the excel- 
lent agreement between the experimentally measured S 
(filled circles) ^'^ and that calculated numerically (solid 
line) for bulk PbTe. In addition, one can see that for any 
nanoinclusion radius, the Seebeck coefficient is always 
increased compared to that of the inclusion-free system. 
In fact, for the smallest nanoinclusion radius considered 
here (1.5nm), the enhancement in S is over 100% at high 
doping. 

It is interesting to consider the impact of Vb on the 
calculated Seebeck coefficient. Figure 3a shows 5 as a 
function of Vb- It is clearly seen from this figure that 
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FIG. 2: Calculated Seebeck coefficient for PbTe with metallic 
nanoinclusions as a function of the doping for several different 
values of the nanoinclusion radius. 



the presence of an extended electrostatic potential leads 
to an increase in S regardless of the sign of Vb ■ {Vb — 
is equivalent to bulk PbTe with nanoinclusions. Nega- 
tive values correspond to the situation of Fig. 1, and 
positive values represent a Schottky barrier). This gen- 
eral behavior can be understood (at least for small |Vb|) 
from the Born approximation which predicts that the in- 
verse scattering time is proportional to the square of the 
perturbation potential. With increase of | Vb | the contri- 
bution to the total inverse relaxation time from inclusion 
scattering increases, leading to an increase of S because 
the energy dependence of r changes from that of T-^uik 
to the more strongly energy dependent r^ . For large val- 
ues of jVel the contribution of island scattering becomes 
dominant and S saturates as seen in Fig. 3a. 

Figure 3a also shows the calculated values of ct as a 
function of the interface potential Vb ■ The conductivity 
decreases as \Vb\ is increased, with a fairly symmetric 
behavior for ±Vb- Combining the results for a and S, 
we obtain the power factor S^a as depicted in Fig. 3b. 
There, one can see that the power factor is increased 
compared to that at Vb = 0, in a range of interface po- 
tentials — O.lSeV < Vb < O.lSeV. The power factor has 
two maxima at some optimal values of Vb because the 
Seebeck coefficient saturates for large | Vb \ while the elec- 
trical conductivity a continues to decrease with increase 

0i\VB\. 

For the optimal interface potential Vb ~ iO.OTeV, we 
find that the power factor is increased by ~ 35%. Impor- 
tantly, the power factor does not decrease substantially 
over a wide range of values of the interface potential. 
Thus, it is possible to take full advantage of the reduc- 
tion in thermal conductivity due to phonon scattering at 
the nanoinclusions, as we will discuss later. 

Instead of optimizing the value of the interface poten- 
tial Vb to achieve the maximum power factor as shown in 
Figure 3b, we can keep Vb fixed (by choosing a specific 
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FIG. 3: Panel (a) shows the calculated Seebeck coefficient and 
conductivity for PbTe as a function of the interface potential 
Vb- Panel (b) shows the resulting power factor. Parameters 
are R = 1.5nm, T = 300K, x = 5%, and n - 



metal for the inclusions) and optimize other parameters, 
for example the inclusion volume fraction or radius. To 
analyze the dependence of the transport coefficients on 
these parameters we can use Eq. ((39)) . As we noted 
above, the function a{E, R) in Eq. ([3D]) varies slowly 
with both E and R, so the inverse relaxation time due 
to electron scattering by inclusions can be approximated 
as T~^{E) « CE^^/^, where the constant C depends on 
Vb, X, and R mostly through the combination (at least 
for small Vb) 



C oc V^x/R. 



(41) 



lir turn, the transport coefficients (S*, a, and Kg) depend 
only on this ratio of parameters. This means that if any 
two parameters out of these three are fixed one can al- 
ways adjust the third parameter (for example, the one 
that can be most easily tuned in the experiment) to max- 
imize the power factor. 

Before closing this section, we remark that we have 
searched for resonant tunneling states in the positive Vb 
regime as a way to enhance S. However, we have not 
found any significant increase in S beyond that already 
discussed. The reason is that, because the potential con- 
tains contributions from several Legendre polynomials /, 
the appearance of a resonant state for one value of I is 
washed out by the nonresonant conditions in the other 
channels. 
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FIG. 4: (color online) Thermoelectric coefficients as a func- 
tion of temperature calculated for n = 2.5 x 10^ ctti^ , 
X = 5%, and R = 1.5 nm. (a) The optimized ZT factor, with 
the optimal values of Vb shown in the inset. (b) The Seebeck 
coefficient, (c) The electrical conductivity. Inset shows opti- 
mized (filled circles) and bulk PbTe (solid line) power factor 
S^a. (d) The thermal conductivity. In all panels, solid cir- 
cles include electron and phonon scattering on the inclusions; 
open circles include electron scattering on the inclusions but 
with bulk PbTe values of Kph', solid lines and filled squares 
are the calculated and measured^^ values for bulk PbTe. 



B. Enhancement of the ZT factor 



While S and a are quantities of fundamental inter- 
est, for applications it is usually ZT = S'^aT/n that is 
most important. The electronic relaxation time calcu- 
lated above leads directly to the electronic thermal con- 
ductivity Kg. Since the total thermal conductivity is the 
sum of electronic and phonon contributions, k — K^+Kph, 
to obtain ZT we also need to calculate Kph- For this pur- 
pose, we adopt a previous method^' ^^ that considered the 
scattering of phonons on nanoinclusions, with the scat- 
tering mechanism for short wavelength phonons being the 
different sound velocities in the host and nanoinclusions. 
This approach has been shown to give excellent agree- 
ment with experiments on nanoscale ErAs inclusions in 



InGaAs^. For T > To {To = ISOiiT is the Debye tem- 
perature of PbTe) Kph can be written as^ 



^ph 



2Tr^v,h^ 



r 

J c 



° Tph {fiuj) d (huj) 



(42) 



where Vg is the speed of sound in PbTe and hiu is the 
phonon energy. The phonon relaxation time Tph is given 

by 



' ph 



• D 



(43) 



where 



cTw^ is the contribution of umklapp 



scattering^ and td is due to scattering by nanoinclu- 
sions. The constant c was determined from Eq. (|42| us- 



ing the experimental value k^J^^ = 2.0W/mK for PbTe 
at T = 300K. For td we used the expression derived in 
Refs. [I6lll4|. In the near geometrical scattering regime 
{qR ^ 1) T£) reads 

T^i - n,Vs{27TR^)[l - sin(2e)/e + sin^{0/C% (44) 

where ^ — qR{va/v'g — \), q is the phonon wave vector, and 
v'g is the speed of sound inside the inclusion. Numerical 
tests show that when the difference in the sound velocities 
is larger than 20% the integrated quantity Kph weakly 
depends on this difference and td can be approximated 
by its geometrical limit value 



niVs{2nR 



3 X 
2r' 



:Ws- 



(45) 



Note that this phonon scattering regime is opposite to 
that on point defects where qR ^ 1. 

Figure 4 shows ZT and its components calculated for 
X = 5%, R = 1.5nin and a doping n = 2.5 x 10^^cm~^, as 
a function of temperature. We discuss this doping first 
because experimental values of ZT and all of its compo- 
nents are readily available for inclusion-free PbTe, and 
can be compared with our calculations; indeed, the calcu- 
lated values for ZT, S, a, and k (solid lines in Fig. 4a-d) 
are in good agreement with experiments^ (filled squares) 
for T < TOOK. (The deviations for T > 700K origi- 
nate in our neglect of the hole contribution to the charge 
and heat transport.) In the presence of nanoinclusions, 
the individual components of ZT deviate from their bulk 
PbTe values at all temperatures shown. For T > AOOK 
the increase of the Seebeck coefficient is compensated by 
decrease of the conductivity, and the optimized power 
factor is close to that of bulk PbTe (see inset in Fig. 4c). 
At such temperatures the small increase of the ZT factor 
due to 'electron-only' scattering by nanoinclusions (open 
circles in Fig 4a) is a result of the decrease of Kg (open 
circles in Fig. 4d). Comparing the ZT shown by filled 
and open circles in Fig. 4a one can conclude that at a 
doping n — 2.5 x lO^^cm"'^ the enhancement of the opti- 
mized ZT is primary due to decrease of Uph , at least for 
T > 400 K. 

To get a more comprehensive understanding of the role 
of nanoinclusions in enhancing the thermoelectric prop- 
erties, we show in Fig. 5 the room-temperature ZT factor 



as a function of the interface potential for two values of 
the doping. In addition, we plot ZT calculated using the 
bulk value of the phonon thermal conductivity k^^"^ (dot- 
ted lines). We first consider the situation of high doping, 
as depicted in panel (a). In the absence of a spatially- 
varying potential [Vb = 0) and without phonon scatter- 
ing on nanoinclusions, ZT is given by the filled circle. 
Turning on the phonon scattering gives a modest 25% 
increase in ZT (the star in the figure). Similarly, one 
can consider the increase in ZT without phonon scatter- 
ing on the nanoinclusions (dotted line); in this case, a 
large increase in ZT of up to 224% is obtained. Thus at 
this doping, electron scattering can give a much larger 
increase in ZT. However, the true advantage of nanoin- 
clusions is realized when both electron and phonon scat- 
tering are included, and the ZT factor can be increased 
by as much as 430%. This increase is much larger than 
simply the sum of the individual electronic and phonon 
contributions. 

The origin of this behavior lies in the non-additive 
effects of electronic and phonon thermal conductivities, 
since ZT depends inversely on their sum. For the large 
doping situation of Fig. 5a we have Kg"'*^ = A.2W / raK , 



bulk 



2.0W/mK, and therefore k! 



bulk 
■ph 



< < 



bulk 



In this 



case, reducing Kph by itself does not lead to an apprecia- 
ble gain in ZT. However, when Ke is also reduced be- 
cause of scattering and becomes comparable to Kph, then 
both work in concert and lead to a large increase in ZT. 
Thus, one can imagine that electron scattering on the 
electrostatic potential serves as an amplification mecha- 
nism to enhance the impact of the reduction in phonon 
thermal conductivity. This mechanism works here be- 
cause at high doping (1) k'^^'' < k^""^ and (2) the power 
factor is maintained or even enhanced in a wide range of 
interface potentials. 

The situation is quite different in the case of low dop- 



ing, where k\ 



bulk 
ph 



> K 



bulk 



as illustrated in Fig. 5b. In this 



case, the electronic thermal conductivity is already quite 
low, K^""= = 0.6W/mK < nlf" = 2.0W/mK, and the 
main impact of nanoinclusions is to decrease the phonon 
thermal conductivity. The maximum increase in ZT is 
107%, with 94% coming from phonons alone. In fact, for 
this low doping, the power factor is always reduced com- 
pared to the inclusion-free system - a signature of this 
effect is the reduction of ZT below that of the inclusion- 
free system for larger values of Vb ■ 

Fig. 6 shows the calculated ZT as a function of tem- 
perature for low (n = 5 X IO^^ctti^'^) and high (n == 
5 X 10"'^^c77i~^) doping levels. Included in the figure are 
the ZT factor calculated with both electron and phonon 
scattering on nanoinclusions (filled circles), that calcu- 
lated with 'electron-only' scattering on nanoinclusions 
and with bulk PbTe values of Kph (open circles), and the 
ZT calculated for inclusion- free bulk PbTe (solid line). 
The corresponding values of Vb that maximize ZT are 
shown in the insets by filled circles. The inset in Fig. 6a 
shows that at n = 5 x lO^^cm"'^ the optimal Vb is very 
small \Vb\ < O.OSey and even vanishes for T > 600K. 
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FIG. 5: Temperature dependence of the optimized ZT fac- 
tor for PbTe. (a) low doping n = 5 x 10^*cm~^ and (b) 
high doping n = 5 x lO^^cm,'^. In both panels filled circles 
denote the optimized ZT calculated with both electron and 
phonon scattering on nanoinclusions, open circles denote ZT 
calculated with electron scattering on nanoinclusions and with 
bulk PbTe values of Kp^, and the solid line is for bulk PbTe. 
The inset in (a) shows the values of Vb that maximize ZT . 
In (b) the bottom inset shows the optimal values of Vb that 
maximize ZT (filled circles) and Vb^ for Pb nanoinclusions. 
The top inset shows the calculated power factor. 



Thus, the electron contribution to enhancement of op- 
timized ZT is negligible (solid line and line with open 
circles almost coincide in Fig. 6a) and the enhancement 
of the optimized ZT is dominated by the reduction in 
Kph due to phonon scattering on the inclusions. This can 
be explained by the fact that for bulk PbTe the Seebeck 
coefficient increases with decrease of the doping concen- 
tration n, and the relative enhancement of the Seebeck 
coefficient from its bulk value due to electron scattering 
on inclusion is smaller at low doping compared to high 
doping (see Fig 2). As a consequence the reduction of 
u at low doping overweights the increase of S"^ and the 
power factor is reduced compared to the inclusion-free 
system, leading to small or vanishing optimal Vb- From 
a practical point of view this result means that in order 
to enhance the ZT factor at low doping levels one needs 
to find a metal that gives little or no interfacial potential. 

For larger doping the electron contribution to enhance- 
ment of ZT becomes important and the optimized Vb in- 
creases. It is seen in Fig. 6b that for n = 5 x lO^^cjn^"^ 
the 'electron-only' contribution to enhancement of the 
optimized ZT makes up over 50% of the enhancement at 
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FIG. 6: Panel (a) shows the calculated Seebeck coefficient and 
conductivity for PbTe as a function of the interface potential 
Vb- Panel (b) shows the resulting power factor. Parameters 
are R = 1.5nm, T = 300K, x = 5%, and n = 2.5 x 10^^ cm"^. 



T < 600K and the optimized Vb is as large as Q.2eV . 
Moreover, at large doping, k^""^ > k^^'*^, and the reduc- 
tion of Ke due to scattering on inclusions amplifies the 
effect of the reduced Kpu- The upper inset in Fig. 6b 
shows that the power factor aS"^ is enhanced only for 
T < QOOK. The reduction of the power factor relative 
to the inclusion-free system at T > 600K is due to the 



fact that at high temperature the Seebeck coefficient of 
bulk PbTe is already large (see Fig. 4b), the relative in- 
crease in S induced by electron scattering on inclusions 
becomes smaller at increased temperature, so the reduc- 
tion of a overweights the increase of S*^. Nevertheless, 
the 'electron-only' contribution results in enhancement 
of ZT (open circles in Fig. 6b) even at higher tempera- 
tures due to reduction of Kg. 

Figure 6b shows by open squares the ZT factor 
for PbTe with Pb inclusions (assuming Vg'' — Ep = 
— 0.35ey) for parameters n = 5 x lO^^cm"^, R = l.bnm, 
and X = 5%. For this set of parameters the interface po- 
tential Vg^ is close to the optimal one (see inset in Fig 
6b) in a wide range of temperatures, so the ZT factor for 
the system with Pb inclusions is very close to the opti- 
mal ZT . The enhancement of ZT due to Pb inclusions 
is on the order of 400% at room temperature and 50% at 
T = 900 K, where it reaches a value as high as 1.5. 
V. CONCLUSION 

In conclusion, we developed a theory that allows the 
calculation of the ZT factor and its components for a 
system of a semiconductor host material with spherical 
metallic nanoinclusions. The enhancement of the See- 
beck coefficient can be explained by a strong energy de- 
pendence of electron scattering on the band-bending at 
the interface between metallic inclusions and the semi- 
conductor host. The electronic contribution to enhance- 
ment of ZT is important for high doping, while at low 
doping the enhancement of ZT is dominated by the re- 
duction in the phonon thermal conductivity. The theory 
can be used to choose the optimal parameters for the 
metal nanoinclusions (interface potential, inclusion vol- 
ume fraction or radius) in order to maximize ZT. 
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